Lifetime measurements of electron spins in Si/SiGe quantum dots 
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We have observed the Zeeman-split excited state of a spin-Y2 multi-electron Si/SiGe depletion 
quantum dot and measured its spin relaxation time Ti in magnetic fields up to 2 T. Using a new 
step-and-reach technique, we have experimentally verified the (/-value of 2.0 ±0.1 for the observed 
Zeeman doublet. We have also measured T\ of single- and multi-electron spins in InGaAs quantum 
dots. The lifetimes of the Si/SiGe system are appreciably longer than those for InGaAs dots for 
comparable magnetic field strengths, but both approach one second at sufficiently low fields (< 1 T 
for Si, and < 0.2 T for InGaAs). 



Introduction Gate-defined quantum dots (QD) in 
the Si/SiGe material system have been touted as one 
of the more promising candidates for spin-based quan- 
tum computation, primarily because of the long deco- 
herence time T 2 expected for electron spins in Si pQ. 
Because T\ establishes a natural upper bound on T2, a 
measurement of the spin relaxation time T\ for isolated 
electrons in any Si-based system has been an impor- 
tant and actively-pursued objective [21 E]. Although 
there has been significant progress in fabricating and 
testing few-electron Si/SiGe dots j3], there has been, 
up until now, no direct measurement showing the an- 
ticipated long spin-relaxation lifetimes. Indeed, there 
has been no direct measurement of a Zeeman splitting 
for either single- or few-electron states in any Si-based 
quantum dot. 

In this Letter, we report the first direct experimental 
confirmation of a Zeeman spin excited state in a few- 
electron Si/SiGe quantum dot, and a measurement of 
its lifetime as a function of magnetic field. We also 
report comparable data for InGaAs dots, used as a 
testbed for our measurement techniques. These dots 
could be emptied, thus allowing us to measure single- 
electron spins. 

Dot Particulars Figure [l] shows the electrode lay- 
out and epitaxial structure for our Si/Sio.7Geo.3 dot. 
The InGaAs/InAlAs/InP depletion dot had an almost 
identical electrode geometry, but a different epitax- 
ial structure which incorporated two quantum wells, 
not one, so that the wafer could also be used for 
accumulation-mode devices (described elsewhere jj]). 
In our depletion-mode devices, the upper well in the 
InGaAs structure played no role. 

Measurement The dot was electrostatically 
formed using conventional techniques [5]. We found, 
empirically, that the lowest number of electrons on the 
dot, N m i n , can be reached when the T-R channel is 
pinched off and the T-B channel kept open. However, 
as the electrode voltages are lowered to drive additional 
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Figure 1: (a) A SEM image of the metal electrodes on the 
device surface. The dashed circle is the estimated location 
of the charge disc formed by the dot electrons. The ar- 
row shows the tunneling channel used in all measurements, 
(b) The CVD-grown Si/Sio.7Geo.3 heterostructure with a 
2DEG mobility of 10 000-15 000 cm 2 /Vs and a charge den- 
sity of 4.6-5.3 x 10 11 cm~ 2 (both measured at 4.2 K). 



electrons off the dot, the T-B channel eventually be- 
comes pinched off also, preventing any further reduc- 
tion in N. Although we routinely reached N — in 
InGaAs dots before this happened, the lowest value of 
N obtained for Si/SiGe dots was likely 7, and this only 
after several painstaking redesigns of the electrodes in 
which the gaps for the T-B and T-R channels were 
progressively widened. 

When iV m i n was reached, T\ was measured by re- 
peatedly applying a three-step bias sequence to one 
of the electrodes (a technique first used in [6]). The 
response of the charge-monitoring currrent, Iqpc, to 
the bias changes is twofold: (i) a drop in the current 
by a fixed amount when an electron is (spontaneously) 
added to the dot, and an increase when it is removed, 
as shown by the dashed lines in Figure § and (tf) a 
replication of the bias sequence due to unwanted ca- 
pacitive coupling of the driven electrode to the R-Q 
channel. 

The key response is the rectangular pulse formed 
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Figure 2: The Delft three-step sequence used for the Ti 
measurement [6]. The first step loads an electron into the 
excited state with probability po, or into the ground state 
with probability (1 — po), and then waits. If the electron 
is in the excited state, and survives the wait-time without 
decaying to the ground state, the second step removes it 
from the dot and reloads another from the reservoir into 
the ground state. The third step flushes the ground state 
electron off the dot, so that the cycle can be repeated. 
The dashed lines are an artistic rendition of instantaneous 
values of 7qpc that occur during a typical cycle, the solid 
curve is an actual time-average of 10 000 of these events. 



during the read interval, which only forms when a 
loaded excited-state electron does not decay to the 
ground state during the wait-time. A time-average of 
this rectangular pulse gives a "spin bump" having the 
analytic form 

h(t) ocpie-t^ - e-'/^/^/n - 1) 

where t\ (t 2 ) is the tunneling time off of (onto) the 
dot, and p is the fraction of electrons that load into 
the excited state during the inject pulse and survive 
the wait-time without decaying to the ground state. 
For equal tunneling times this expression reduces to 
h{t)ocp{t/T)e- t / T . 

With increasing wait-time i wa it, the fraction of in- 
stances when an unload-reload pulse is formed dur- 
ing the read interval will decay exponentially as p = 



ait/Tl 



(assuming that the characteristic time 
Ti ^> r). The spin bump height /i max , plotted versus 
twait, will show the same exponential decay, as demon- 
strated by the data in Fig. [3] 

There are two advantages of using time-averaging 
rather than the pulse-counting of [6]: (i) the high 
signal-to-noise ratio required to reliably detect indi- 
vidual pulses is not needed, and (it) the false counts 
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Figure 3: Spin bump height versus wait-time for our 
Si/SiGe dot in a magnetic field of 1.5 T. The solid curve, an 
exponential fit to the data (circles), gives a Ti of 164 ms. 



due to random telegraph signals (RTS) that will un- 
avoidably occur when the read level is near the lower 
Zeeman level contribute trivially, because the average 
of RTS-induced pulses is a straight line. What is lost 
by using the averaging approach is a direct measure of 
the branching ratio, p . However, p is not needed for 
a determination of T\. 

Zeeman splitting The spin-bump appears when, 
and only when, the amplitude of the inject-and-wait 
pulse is sufficient to reach the first excited state. This 
suggests a step-and-reach technique for measuring the 
Zeeman splitting. One places the read level a known 
amount above the ground state (the step process), and 
then progressively reaches up with the inject pulse un- 
til a bump appears. The exact point of appearance 
of the bump is determined experimentally by measur- 
ing the bump height for several values of step-plus- 
reach amplitude just beyond the first appearance of 
the bump, and then extrapolating backward to de- 
termine the zero- intercept. To convert the combined 
step and reach amplitudes to an actual energy, we 
used a temperature sweep and the known relation be- 
tween the width of a Coulomb Blockade peak and the 
dot's temperature [7] to determine both the lever arm, 
a, and the effective base temperature of the system, 
12.55 Vpi un ger/eV and 48 mK, respectively. 

The Zeeman splitting was measured for several val- 
ues of the magnetic field using this technique. Two 
sets of measurements were made, one with the mag- 
netic field increasing with each successive point, the 
other with it decreasing, in order to look for hysteretic 
effects (which we did not see). These data are plot- 
ted in Fig. [4] One sees immediately that there is a 
small but non-negligible deviation from a straight line 
at higher values of B, which we suspect may be caused 
by a small component of the magnetic field normal to 
the 2DEG. We have fit the 9 measured points to the 
somewhat arbitrary function Ez = a + bB + cB n , using 
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Figure 4: The Zeeman splitting as a function of magnetic 
field. The solid curve is a fit with Ez — a + bB + cB 3 , the 
dashed curve is the linear portion of that expression. 

values of n ranging from 2 to infinity, and have found 
that these fits give a zero-field value of g that ranges 
from 2.1 to 1.9. Although the g-value that results from 
the fit having the lowest variance and proper symme- 
try (n = 3) is 1.99, all fits appear equally acceptable 
to the eye, making it difficult to choose one over an- 
other. Our estimate of g is thus 2.0 ± 0.1, in agree- 
ment with theoretical expectations and measurements 
on impurity-bound electrons in Si [8 . 

Multi-electron dots Our measurements of T± for 
Si/SiGe are for an iV-electron dot that we believe had 
7 electrons (5 to 7 electrons if based exclusively on the 
energy spacings between Coulomb Blockade peaks, or, 
more definitively, 7 electrons when further constrained 
by the filling sequence of spin states as discussed be- 
low). Measurements of T\ at lower values of N were 
not possible because the dot "closed". 

To see what would happen at values of N > 1 in dots 
that can reach N — 0, we looked for spin bumps on 
our InGaAs dot, starting with N = 1, and progressing 
all the way up to N — 9. We found a spin bump 
at N = 1, 3, 5, and 7, but none for 9, and none for 
N = 2, 4, 6, and 8. These sightings are consistent with 
Pauli filling of an asymmetric dot |9j, in which the spin 
values would be 1/2, 0, 1/2, ... for N = 1, 2, 3, 4. . . 
electrons. The measured value of Ti for N — 7 was 
12.8 ms, which is shorter than the 41.5 ms measured 
for N = 1, but not dramatically so. 

For Si/SiGe dots, however, which have an extra de- 
generacy due to valleys, there is reason to believe (e. 
g., see [TO]) that the filling sequence for an asymmet- 
ric dot should be 1/2, 1, 1/2, 0, 1 /2, 1, 1/2, . . . when the 
Zeeman energy exceeds the valley splitting. For our 
lowest N-value, the measured Zeeman splitting and g- 
value of 2.0 imply a spin of Y 2 ( or greater). If the 
spin is Y2, the next higher value of N should have an 
s = or s = 1 ground state. Measurements of the 
N + 1 dot showed that the ground state had s — 0, 



and the excited state 8 = 1, with a surprisingly-small 
singlet-triplet splitting of approximately 36 fj,eV jll] . 
Thus we contend that our measured values of T\ are 
for an unpaired spin of Y2, and are representative of 
(although probably somewhat shorter than) the value 
that would be obtained for a single-electron spin- 1/2 
dot. 

Theory At reasonably strong spatial confinement, 
low temperatures, and small-to-moderate magnetic 
fields, the single-phonon admixture mechanism dom- 
inates the intrinsic spin flip (Ti relaxation process) of 
zero-dimensional conduction electrons in III-V mate- 
rials |12j . It can be understood as follows. Spin-orbit 
(SO) coupling admixes higher orbital states with oppo- 
site spin projections into the eigenstates of the lowest 
spin doublet, thus allowing, in principle, intra-doublet 
transitions (i. e., spin flip) to be induced by strictly 
spin- independent interactions with the environment, 
e. g., with lattice vibrations. Exact cancellation in 
the transition matrix element, enforced by the time- 
reversal symmetry of the total electron Hamiltonian 
(including the SO admixing terms), should be broken 
by a magnetic field to yield a finite spin relaxation rate. 

At low magnetic fields the Zeeman splitting of the 
SO-mixed doublet (and, thus, the required energy 
transfer to the lattice phonon) is small, so that, when 
present, piezoelectric electron-phonon coupling domi- 
nates, resulting, at temperatures T <C gfisB /kg, in a 
well-known |12l [T3l H"4"] functional dependence 

T-p 1 - P V 2 B 5 /El 

obtained in the (valild for our case) approximation 
of long- wavelength resonant phonons. Eq is the spa- 
tial quantization energy. SO coupling is assumed to 
be dominated by linear-in-fc SO terms in the electron 
dispersion; their strength is quantified by a material- 
and structure-specific constant 77, which is, in general, 
anisotropic. In non-piezoelectric materials, but also 
in piezoelectric materials at larger magnetic fields, en- 
ergy transfer to phonons is facilitated by a deformation 
potential coupling resulting (e. g., [13]) in 

T-J = DrfB 7 /E* 

at low temperatures. At temperatures T 3> g^sB/kB, 
the spin-flip rate is given by ksT/ gfisB. 

Moderate electric fields < 2 mV/nm are typi- 
cally present in our dots. In InGaAs structures, 
a bulk-inversion-asymmetry-induced (BIA or Dres- 
selhaus) SO term with j3 ~ 1 600 /ieV nm dom- 
inates the structure- induced (SIA or Rashba) term 
(< 500 /ieV nm), making r\ ss (3 a field-independent 
and isotropic constant, so that no in-plane anisotropy 
is to be expected for T% in a circular dot. P and D are 
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main crystallographic axes As formulated, these 
mechanisms are not directly related to the heterointer- 
face properties. They should exceed T^ 1 at extremely 
small magnetic fields, but the estimated times are too 
slow to be observed in our experiments. 
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Figure 5: Inverse spin lifetimes as a function of magnetic 
field. Circles are for the Si depletion dot, squares (dia- 
monds) for the InGaAs depletion (accummulation) dot. B 
is always parallel to the 2DEG (along the [010] axis). The 
straight lines are theoretical predictions for the two mate- 
rial systems. 

material-, but not structure-specific parameters (apart 
from a weak g factor dependence). When expressing 
Tj~ in inverse seconds, rj in fieV nm, the in-plane B in 
Tesla, and E$ in meV, the InGaAs-specific constants 
P and D are found to have numeric values of 7.6 and 
0.26, respectively |15j . These low-temperature asymp- 
totics of T^ 1 and are shown in Fig. [H] by dashed 
and dash-dotted lines, respectively, for the relevant 
Eq = 10 meV. The dotted line is the high tempera- 
ture £? 4 -asymptote of T^> calculated for T = 1 K, Our 
measurement temperatures are always much lower. 

The Si crystal symmetry allows neither piezoelectric- 
ity nor a BIA SO term. SIA is dwarfed by an especially 
weak bulk SO coupling for A valleys (estimated to be 
only 3-4 /^eV nm in comparable structures [16]). It was 
suggested in [T7] that linear-in-fc SO terms in (001) Si 
heterostructures are dominated by contributions of a 
low-symmetry heterointerface. For our typical electric 
fields, we estimate, with the help of Fig. 3 in [TJ], that 
•q has a maximum value of 20-30 /ieV nm for a per- 
fect interface, but a much lower value in the presence 
of interface imperfections. In Fig. [H] T^j in Si (with 
D = 4.3 x 10~ 3 ) is shown by solid lines for (the nu- 
merically simulated) Eq = 2 meV. 

Other spin-lattice relaxation mechanisms were pro- 
posed for Si dots, in particular: (i) the modulation of 
the hyperfine coupling by phonon deformation, which 
has a B 3 field dependence |12llT8] . and (it) the phonon- 
induced modulation of the bulk electron g-factor, with 
a .^-dependence and a strong sensitivity to the orien- 
tation of the in-plane magnetic field with respect to the 
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